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Abstract 

Within the minimal supersymmetric standard model (MSSM), the large tan/5 regime 
• ^ ! can lead to important modifications in the pattern of CP- violating sources contributing 
^ I to low energy electric dipole moments (EDMs). In particular, four-fermion CP- violating 
' interactions induced by Higgs exchange should be accounted for alongside the constituent 
EDMs of quarks and electrons. To this end, we present a comprehensive analysis of three 
low energy EDM observables - namely the EDMs of thallium, mercury and the neutron - 
at large tan /3, in terms of one- and two-loop contributions to the constituent EDMs and 
four-fermion interactions. We concentrate on the constrained MSSM as well as the MSSM 
with non-universal Higgs masses, and include the CP-violating phases of fi and A. Our 
results indicate that the atomic EDMs receive significant corrections from four-fermion 
operators, especially when Im(y4) is the only CP-violating source, whereas the neutron 
EDM remains relatively insensitive to these effects. As a consequence, in a large portion 
of the parameter space, one cannot infer a separate bound on the electron EDM via the 
experimental constraint on the thallium EDM. Furthermore, we find that the electron EDM 
can be greatly reduced due to the destructive interference of one- and two-loop contributions 
with the latter being dominated by virtual staus. 



1 Introduction 



Electric dipole moments (EDMs) of the neutron [1] and heavy atoms and molecules [2, 3, 
4, 5, 6, 7, 8] are primary observables in testing for flavor-neutral CP violation. The high 
degree of precision with which various experiments put limits on possible EDMs translates 
into stringent constraints on a variety of extensions of the Standard Model at and above the 
electroweak scale (see, e.g. [9]). The experiments that currently champion the best bounds 
on CP-violating parameters are the atomic EDMs of thallium and mercury and that of the 
neutron: 

|(iTi| < 9 X 10"^^ecm 

I^Hgl < 2xl0-2^ecm (1.1) 
\dn\ < 6xl0"^^ecm. 

In this paper we address the theory of electric dipole moments in supersymmetric 
(SUSY) models with a large ratio of Higgs vacuum expectation values, or tan/3. Mod- 
els with large tan f3 have recently received signiflcant attention, stimulated in part by the 
final LEP results that impose significant constraints on the Higgs sector of the MSSM and 
imply a relatively large tan/3, tan/5 > 5 [10]. Another motivation for large tan/3 stems 
from unified theories, which allow for the unification of top and bottom Yukawa couplings 
when tan/3 ~ 50 [11]. 

Supersymmetric models with the minimal field content (MSSM) allow for the presence 
of several CP-violating phases even in the most restrictive ansatz of flavor universality in 
the squark and slepton sectors. The null experimental EDM results pose a serious problem 
for the MSSM with superpartner masses around the electroweak scale. Indeed, a typical 
CP- violating SUSY phase of order one combined with C(100 GeV — 1 TeV) masses for 
superpartners would violate experimental constraints by up to three orders of magnitude 
[12]. This generally requires the CP- violating phases to be very small, unless there are 
cancellations among different contributions [13, 14, 15] but these are largely disfavored by 
the mercury EDM constraint [16, 17, 18]. In string models, there are further complications 
since large EDMs are generally induced even if the SUSY breaking dynamics conserve CP (in 
the sense that the SUSY breaking F-terms are all real) [19]. Suppression of these unwanted 
effects requires rather special circumstances, e.g. dilaton domination with an approximate 
axial symmetry [20]. Other options to avoid overproduction of EDMs in SUSY models 
include heavy sfermions [2 1] , the presence of extra dimensions [22] , additional symmetries 
[23, 24], etc. Yet, despite the sizeable literature on EDMs in supersymmetric models, EDM 
constraints at large tan /3 remain poorly explored. 

With this situation in mind, it is natural to ask what changes to the EDM predictions 
will occur at large tan/3? One simple observation is that the EDMs of the down quarks 
and the electrons, as induced by SUSY CP-violating phases, tend to grow because of their 
dependence on the down-type Yukawa couplings which translates into a linear enhancement 
in tan /3. The second and in some sense more profound change was for a long time overlooked 
in the SUSY-EDM literature. At large tan /3, the observable EDMs of neutrons and heavy 
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atoms receive contributions not only from the EDMs of the constituent particles, such as 
electrons and quarks, but also from CP-odd four-fermion operators. The latter operators, 
since they are induced by Higgs exchange, receive an even more significant enhancement by 
several powers of tan /3. The necessity of including these operators in the EDM analysis was 
shown by Barr [25] for the case of a two Higgs doublet model with spontaneous breaking of 
CP, and recently by Lebedev and Pospelov [26] for the MSSM with exphcit CP violation. 

In this paper, we present a thorough analysis of EDMs in the large tan /3 regime. We 
analyze the predictions for the three obscrvables that champion the best constraints on 
the CP-odd sector of the theory: the atomic EDMs of thallium and mercury [2, 3], and 
the neutron EDM [1]. The restrictive assumption of Ref. [26] on the parametrically large 
superpartner mass scale will be lifted in this paper, and the effects of the four-fermion op- 
erators are analyzed alongside the usual electron and quark EDM and color EDM (CEDM) 
contributions. We first consider a framework in which the scalar and gaugino masses, mo 
and mi/2, as well as the tri-linear terms Aq, are unified at the GUT scale, i.e. the con- 
strained MSSM (CMSSM) [27, 28]. In the CMSSM, the Higgs mixing mass, fi, and the 
pseudoscalar Higgs mass, m^, are determined by the GUT scale parameters, mi/2,mo, Aq, 
and tan (3 using the radiative electroweak symmetry conditions. The sign of // is also free 
in the CP-conserving version of the CMSSM. We then repeat the same analysis for the 
case when scalar universality is relaxed in the Higgs sector (NUHM) [29, 30]. In this case, 
the Higgs soft masses mi and m2 can be chosen independently of the squark and slepton 
masses, which remain unified at the GUT scale. Alternatively, one may choose /x and m^ 
as free parameters and use the RGEs to determine mi and m2 at the GUT scale. The 
positivity of mi and m| at the GUT scale may impose a restriction on the parameter space 
in this case. 

The universality of the gaugino masses at the GUT scale limits the number of physical 
SUSY CP- violating phases to two^, which in the real gaugino mass basis can be identified 
with the phases of the /i and A parameters. We perform the analysis separately for these 
two phases, as the scahng of the EDMs with tan P is quite different in each case. 

The paper is organized as follows. In the next section we review the structure of the 
effective CP-odd Lagrangian at 1 GeV. We argue that at large tan /3 it should be modified 

by the presence of four-fermion operators. We then estimate the SUSY parameter regimes 
where these contributions become important sources for EDMs. Section 3 calculates the 
supersymmetric threshold corrections that induce CP violation via Higgs exchange. Section 
4 addresses the modifications of the thallium, mercury and neutron EDMs due to the 
presence of four-fermion operators and provides general formulae for these observables. In 
section 5 we perform a numerical analysis of EDMs at tan/3 — 50 in the CMSSM and 
the NUHM. Section 6 concludes the paper by summarizing the modifications to the EDM 
bounds at large tan f3 and discussing the importance of different contributions to the EDMs. 

^This is a consequence of the V{1)r and Peccei-Quinn symmetries which allow us to eliminate two 
CP- violating phases (those of Bij, and gaugino masses). In the flavor-nonuniversal case, some of the phases 
can also be rotated away by redefining the quark supcrficlds, such that the physical SUSY CP-violating 
phases are reparametrization invariant combinations of the SUSY and SM flavor structures [31]. 
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Some additional details on relevant hadronic calculations are given in appendices A and B. 



2 The structure of the low energy Lagrangian 

We begin by recalling the standard structure of the effective CP-odd Lagrangian normahzed 
at 1 GeV, which is the lowest scale allowing for a perturbative quark-gluon description. It 
includes the theta term, the Weinberg operator, the EDMs of quarks and electrons, and 
the color EDMs of quarks: 

-- dii^,{Fa)^,^-'- di^MGa)^5^ + --- (2.2) 

i=e,u,d,s i=u,d,s 

For a given pattern of supersymmetry breaking dj, dj, and w can be calculated explicitly, via 
one or two-loop SUSY diagrams, and evolved down to 1 GcV using perturbative QCD. On 
the other hand, only the additive radiative corrections to the theta term can be calculated, 
leaving ^ as a free parameter subject to an arbitrary initial condition. In supersymmetric 
models with CP-violating phases larger than ~ 10~^ in the flavor-diagonal sector the only 
reasonable strategy to avoid the strong CP problem is to postulate the existence of the 
Peccei-Quinn relaxation mechanism [32], which removes the theta term from the list of 
contributing operators. We will therefore adopt this strategy and discard 9. When tan/3 
is low, possible four fermion CP-odd operators (e.g. {qq){qij5q) and alike) can be safely 
neglected. Indeed, these operators involve a double flip of helicity associated with an 
additional {mg/vEwY suppression, and this makes these operators effectively of dimension 
8, and thus totally negligible. 

The dependence of the experimentally measured EDMs of neutrons, and paramagnetic 
and diamagnetic atoms, on the Wilson coefficients in (2.2), as calculated at low tan/3, can 
be symbolically presented in the following way: 

dTl = dTl{de) 

dn — dn{0,di,di,w) at low tan/3. (2.3) 

dug = dng{0,di,di,w) 

The dependence of the EDMs of paramagnetic atoms and molecules (such as thallium) 
on the electron EDM alone makes them the most valuable observables for SUSY models, 
as they allow us to link the experimental results directly to a calculable combination of 
CP-violating phases, and slepton and gaugino masses, while the uncertainty of atomic 
calculations enters as an overall factor. In fact, due to this property of paramagnetic 
EDMs, it is customary to quote the limit on d^ directly, skipping the actual experimental 
limit on the EDM of the paramagnetic atom. The EDMs of the neutron and mercury are 
far more complicated as a variety of hadronic (and in the case of mercury also nuclear and 
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atomic) matrix elements are needed in order translate the experimental results into limits 
on a specific combination of the Wilson coefficients. 

The structure of the effective Lagrangian becomes even more complicated if tan [3 is 
taken to be a large parameter. Indeed, the operators induced by Higgs exchange that 
involve down- type quarks and charged leptons grow rapidly with tan/3 [26]. For now, we 
parametrize these operators by coefficients Cjj, 

E Q.iAAMn^^j): (2.4) 

i,j=e,d,s,b 

where i, j run over flavor indices and the second index always indicates the fermion flavor 
that enters (2.4) via a pseudoscalar bilinear. These operators represent a subclass of the 
larger set of all possible CP-odd flavor-conserving four-fermion operators considered in 
[33, 34, 35] . Note that the 6-quark is heavy compared to the scale of 1 GeV and thus can be 
integrated out, producing new dimension 7 operators ■0j'0j(G^^(7'''^'") and '4^ii^^'4>i{G'^^^,G'^^'"') 
and leading to new higher loop contributions to (ij, di and w. With this procedure in mind, 
we prefer to keep the coefficients C^i and Cjf, explicitly so that integrating out the 6-quarks 
is interpreted as part of taking the corresponding matrix element. 

When the contributions from Cjj are not negligible, the relation between the observables 
and the SUSY parameters encoded in the Wilson coefficients of Eqs. (2.2) and (2.4) becomes 
considerably more involved: 

(ixi = dTl{de,Cde,Cse,Cbe) 

dn = dn{9,di,di,w,Cq^q.J at large tan/5. (2.5) 
c^Hg = d}ig{9,di,di,w,Cij). 

Paramagnetic EDMs lose their one-to-one relation to d^, as they also become dependent on 
operators that involve quarks and gluons, and as a consequence are subject to considerable 
hadronic uncertainties. Hence, at large tan/3, a separate bound on d^ simply may not 
exist, and the bound on the EDM of the paramagnetic atom (molecule) as a whole should 
be used instead. The mercury EDM receives a number of additional contributions from 
semileptonic and pure hadronic four-fermion operators. As a consequence, {du — dd) may 
no longer be the dominant contribution to d^g in contrast to the situation at low tan/3 
[16]. Thus, the large tan/3 regime is qualitatively different from tan/3 ~ and requires 
a separate dedicated analysis. 

To conclude this section, we will determine a combination of SUSY parameters that 
regulates which effect is more important: one-loop EDMs or four-fermion operators. The 
required interpolating parameter can be deduced from the asymptotics of the one-loop 
EDMs and four-fermion operators in the limit of a generic heavy superpartner mass scale 
M. For simplicity, and in this section only, we will not differentiate among different super- 
partners and assume the same mass for squarks, sleptons, gauginos, etc. 

Let us first consider a non-zero /x-phase (relative to the gaugino mass) . While the one- 
loop EDMs generally scale as di/mj ~ //mAtan/3sin^^/M^, the leading order contributions 
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to Cij scale as C ~ //mA(tan/3)^sin^^/(M^m^) [26], where is the CP-violating phase, 
m\ is a gaugino mass, and M — max(mo, ^1/2) stands for a sfermion or gaugino mass, 
whichever is larger. Thus, we can introduce the following combination of SUSY parameters 
that interpolates between the two regimes: 

M 

i = tan/3. (2.6) 

rriA 

In the limit of large M, with uia fixed at the electroweak scale, all EDM-like observables 
are dominated by the four-fermion operators [26]. When M is kept finite, there is a certain 
critical value for ^ below which the usual one-loop EDMs dominate the observables, and 
above which the four-fermion operators are more important. 

Using the results of Refs. [26, 16], we can estimate the critical value for ^ in the hmit 
when all superpartner masses arc approximated by some heavy scale mg ~ '021/2 ~ M ^ 
Vew ~ 246 GeV, while mA is kept fixed near the electroweak scale. If all CP violation is 
induced by 6*^, the relative importance of the two sources is determined by the ratio 



10" 



M 

tan p 

rriA 



= (2.7) 



0?Tl(4) 

which suggests a critical value for ^, 

^e~ 300 for e^^Q. (2.8) 

From this estimate it becomes clear that when tua — M the effects of the four-fermion 
operators are subleading but not negligible, reaching 20% of the dg contribution at maximal 
tan/?. On the other hand, even a mild hierarchy between the mass scales, M ^ 5m^, may 
lead to a contribution to atomic EDMs comparable to d^, or even to dominance of the 
four-fermion contributions. 

In the case of the other CP-violating phase 6a, the relative phase of the A parameter 
and the gaugino mass, we estimate that the transition to dominance of the four-fermion 
operators takes place at 

^ctan^/2^ - 500 9aj^0. (2.9) 

Here, the additional power of tan/5 arises from the fact that the one-loop induced de{6A) 
does not grow with tan f3. Thus in this case, even for tua ~ M, the contribution of the 
four-fermion operators may become comparable to d^. 

Both conditions (2.8) and (2.9) can be realized in a generic SUSY spectrum, but how re- 
ahstic are they for the constrained MSSM and the MSSM with non-universal Higgs masses? 
We note that the simple comparison of de and Cij made in this section may be subject to 
significant modifications, as the masses of sleptons and squarks as well as other super- 
partners are affected by the renormalization group evolution from the GUT scale to the 
electroweak scale. Therefore, we now turn to a numerical comparison of all contributions 
to the EDMs for different CP-violating phases and different patterns of SUSY breaking at 
large tan /3, as detailed in the next four sections. 
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3 CP-odd four-fermion operators in the MSSM 



We now turn to specific SUSY mechanisms that generate Cij. As shown in Ref. [26], Cij are 
generated primarily through the exchange of the pseudoscalar and scalar Higgs particles A 
and H with CP violation entering through vertex corrections. This dominates the one-loop 
induced Cij due to large 0{tan^ (3) enhancements. Other contributions, such as CP-odd 
higher-loop effects, SUSY box diagrams, etc., can have an even stronger dependence on 
tan f3. However, these effects are further suppressed either by additional loop factors or by 
additional powers of small Yukawa couplings and in this sense are subleading. Nevertheless, 
A — H mixing may be numerically important [26], and therefore its contribution has to be 
included in the EDM analysis. 

3.1 Higgs-mediated CP-odd four-fermion operators 

The tan/? enhancement of certain loop corrections in the down- type quark and charged 
lepton sectors is a well studied phenomenon [36]. Supersymmetric threshold corrections, 
such as those in Fig. 1, generate a new Yukawa interaction which is absent in the limit of 
exact supersymmetry. As a result, below the superpartner scale, the Yukawa interactions 
in the D-quark and charged lepton sectors are given by 

-Cy = {Y^'^ +Y},'^)H,DiDL + yD hId^Dl 

+ + Yi^)HrEiEL + yE hIe^El + h.c. , (3.10) 

where Y^\ are the tree level Yukawa couplings and yD,E and Yj^\ are the one-loop-induced 
couplings. Thus, instead of being the Yukawa sector of a type II two-Higgs doublet model, 
Eq. (3.10) represents a generic two Higgs doublet model where Higgs exchange may lead 
to violation of CP. One should keep in mind that Eq. (3.10) is vahd in the hmit of heavy 
superpartners. If the scale of SUSY breaking is close to the electroweak scale one should 
also include additional radiatively generated operators, e.g. {H\H2YH\D^j^Dl. 

Due to fiavor universality in the soft-breaking sector, flavor-changing effects can be 
ignored in EDM calculations, and a simple linear relation between and yD,E holds for 
every D and E flavor: 

yD = JDY^^\ yE = JEYi\ (3.11) 

The calculation of the loop functions Jd and Je is straightforward but tedious, as a large 
number of SUSY diagrams contribute. Nevertheless, an analytical result for this set of loop 
diagrams is possible, and we present it in the next subsection. 

SUSY CP violation arises from a mismatch between the phases of the and A parame- 
ters, and mi/2, and leads to complex masses and Yukawa couplings. The quark and lepton 
masses, 

-Cr.^ = {Y§'^ + YS'^){H^)DiDL + yD{H2)DiDL 

+(4°^ + Yi}^){H^)El,EL + yE{H2)El,EL + h.c. , (3.12) 
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(d) 

Figure 1: One- loop SUSY threshold corrections in the down quark and charged lepton 
Yukawa sectors. All four types of diagrams are relevant for the 6-quark, while only (a)-(c) 
contribute in the case of d and s quarks, and the electron. Here Aj denote the gauginos. 
The vacuum insertion of {H2) indicated by a cross is for illustrative purposes only as all 
calculations are performed in the mass eigenstate basis. 




can be made real via a chiral rotation of the D and E fields, 

Dr^c-'^'^^Dr, Er^c-'^'^^^Er, (3.13) 
where the rotation angles are defined by tan/3 and the loop functions Jf. 

= Arg(y/°^ + y/'^ + yi tan /3) ~ Arg(y/°) + y/'^ ) + Arg(l + tan /3) , (3. 14) 

where i = d,s,b,e arc the flavors of interest. Y^^^ contains a loop suppression factor and 
thus can be neglected. However, a large value of tan (3 can compensate the loop suppression 
of yi so that the phases Arg(l + Jj tan l3) can be of order one. The vertices of the A and H 
Higgs bosons and fermions, as derived from the Lagrangian (3.10), also contain a complex 
phase: 

<5f ^ = Arg((r/°) + r/^)) tan/3 + y,) ~ Arg(r/°) + Y^'^). (3.15) 

This equation is written under the assumption that tan /3 is a large parameter. The differ- 
ence between these two phases, 

5, = 5(-)_5f), (3.16) 

constitutes a physical CP- violating phase that enters into the interactions of H and A with 
physical fermions ipi, 

J^CPY = - ^ sin 5i (i^ii/JiA + ^pii'-f^ipiH^ . (3.17) 

i=d,s,b,e V 
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Figure 2: The Higgs- mediated electron-quark interaction Cqe with CP violation at the 
Higgs-quark vertex (a) and at the Higgs-electron vertex (b). 



Note that the phase Arg(F/°^ + Y^^') drops out of the expression for 5i as it should. An 
exchange by the physical Higgses, as in Fig. 2, will then produce CP-odd four-fermion 
interactions with the following coefficients 



CiJ''\i,j = e,d,s,b) 



tan2/3 Yi 



SMySM 
3 



sin((5j 



2m\ 



|1 + Jitan/3||1 + J,-tan/3| 



(3.18) 



Here Y^^ denote the Standard Model values for the Yukawa couplings, Y^^ = ^/2mf/v 

(t'=246 GeV), and the superscript of Cf™"* signifies that this contribution originates from 
a vertex correction. Note also that at leading order in tan (3 an exchange by the (lightest) 
/i-Higgs does not contribute to C^, the difference between and niH can be ignored, and 
the neutral Higgs mixing angle is trivial, cos^ a ~ 1. 

Examining expression (3.18), we notice immediately that 



-<-y(vc) _ 



(3.19) 



and thus vertex corrections do not generate four-fermion operators that involve a single 
flavor such as, for example, dddi'j^d. 

Combining Eqs. (3.14) and (3.18), we observe that the coefficients Cij grow as (9(tan^ /3) 
as long as Jj(j) tan /9 < 1, which is true for almost the entire available domain of tan (3. 
Moreover, the cubic growth of Cij generated at one loop due to explicit CP violation in 
the soft breaking sector is distinct from the 0(tan^ /3)-behavior of C^j in two-Higgs doublet 
models with spontaneous CP violation [25] . 



3.1.1 Beyond the limit of heavy superpartners 

In order to go beyond the approximation of heavy superpartners, one should reformulate 
Eq. (3.18) in terms of the tree-level fermion masses ml = f/°^wi/V2 and the one-loop 
SUSY corrections to the fermion masses Amj, where all operators neglected in Eq. (3.10) 
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Figure 3: Diagrams that generate subleading contributions to Cij. (a) is the radiatively 
generated H — A mixing, and (b)-(c) are typical box diagrams. 



are included^. The one- loop-corrected mass terms can be made real by a phase redefinition 
(3.13) with the phase now given by 



=Arg(l + Am,/mf ) 



(3.20) 



Here we assume for convenience that m-°-* is real. In the large tan (3 regime, the couplings 
of A and H Higgses to fermions can be obtained as a derivative of mf^ + Ami with respect 
to (Hi) = vi at Vi = 0, and thus the phase S^^^ is now given by 



^(H) ^ ^{mf^+Amj) "^ 



(3.21) 



In the large tan/5 regime it is natural to expect that 61 '^61 , so that the physical CP- 
violating phase 6i is given by Arg(l-|-Amj/m-°^). However, in practice it often happens that 
S^^^ cannot be neglected (see Sec. 5.1). Finally, the coefficients 1-|- Jjtan/3 appearing in the 
denominator of Eq. (3.18) should also be replaced by the expression 1 -|- Arrii/mf^ (which 
is independent of mf'^ because Amj is proportional to mf'^). With these modifications, 
Eq.(3.18) takes the form 



CiT\i,j = e,d,s,b) = 



tan^ f3 



Yf^Yf^Um{6,-6j 



2m\ |l + Ami/m, 



(0), 



1 -|- Amj/m 



(0), 



(3.22) 



where the phases 5i are given in the expressions (3.16), (3.20) and (3.21). Eq. (3.22) 
shows that the crucial step in calculating the CP-odd operators is finding the one-loop 
correction to the masses of the rf, s, h quarks and the electron. 

At subleading order in tan/?, there are a large number of diagrams that contribute to 
Cjj, examples of which are given in Fig. 3. The complete calculation of these diagrams 
goes beyond the scope of the present paper. However, as shown in Ref. [26], the CP-odd 

^ A more complete approach would necessitate the use of the Coleman- Weinberg effective potential. For 
our purposes, however, it suffices to truncate the series of non-renormalizable operators. 
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mixing of H — A Higgs bosons may sometimes be numerically important and comparable 
to (3.22). However, since the H — A mixing term contains either a loop factor which is 
not compensated by large tan/5 or a suppression by {Y^^i^Y ia.n (3 , H — A mixing may 
only become significant with the help of an additional mass hierarchy. In particular, if the 
following combination is large [26, 37], 

» 1, (3.23) 



H — A mixing may provide a sizeable contribution to the EDMs. The box diagrams scale 
difi^erently and lose their importance when the superpartners are heavy, itlq ~ "021/2 ^ v. 
Since this is one of the most interesting domains of the parameter space, we include the 
effects ol H — A mixing and neglect the box diagrams. 

The coefficients Cij induced by diagrams 3a can be calculated with the use of the 
following formula: 



(Ag)tan ^/3 ySMy.SM 
~ 2m1 



Cl ' ~ ^ ' - ^ , (3.24) 



where i,j = d, s, b, e and {AH) stands for the one-loop-generated amplitude in the effective 
Higgs Lagrangian, 

JC-h,a = -\m\A' - ^mlH' - (AH) AH + ■■■, (3.25) 

where the (AH) mixing parameter has the dimension of mass squared. Eq. (3.24) also uses 
m% ~ m\ 3> (AH) and tan /? 3> 1. The effect ofA—h mixing, where /i ~ Re (i72 — (-^^2)), is 
negligible under these conditions. In this approximation, it is not necessary to re-diagonalize 
the 3x3 Higgs mass matrix. It suffices to simply account for (AH) by a direct insertion 
on the Higgs line in Fig. 3a. Note that C\f^^ are flavor-symmetric, 

Cjf-) = C^'\ (3.26) 
and, unlike the vertex corrections, can generate a single flavor four-fermion operator. 



3.2 SUSY threshold corrections in the mass sector and scalar- 
pseudoscalar mixing 

In this subsection we provide a complete set of threshold corrections Am^, at leading order 
in tan/3, to the quark and lepton masses (ignoring flavor mixing). Let us first briefly 
summarize our notation which, unless noted below, follows that of Ref . [38] . The chargino 
and neutralino mass matrices are diagonalized by 

C/*M^+y^ = diag(m^+,m^+) , 

N*M^^N^ = diag( m^o, m^o) , (3.27) 
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where C/, V, and N are unitary matrices. In addition, the sfermion mass matrix is diago- 
nahzed by the matrix 



Cf SfC^ 



(3.28) 



that rotates the basis (/l, /r) to the mass eigenstate basis (/i, /2). Here c/ = cos6f, 
Sf = sinOf, and ^/ is the sfermion mixing angle (see [38] for the exphcit form of the mass 
matrices) . In our convention the CP- violating phase appearing in the sfermion mixing mass 
is given by^ 

cl>f^Arg[-{A} + piRf)] , (3.29) 

where Rf — tan P for — —1/2 and Rf — cot /3 for 73 = 1/2. 

As argued in the previous section, a large value of tan/3 can compensate the loop 
suppression of Ami such that the consequent CP- violating phase is of order one. Below we 
list the relevant supersymmetric contributions to Am^ which grow linearly with tan /3. We 
omit the superscript in mf^ for brevity. 

(i) Corrections to the down-type quEirk masses 

For all three down-type flavors, d, s and b (where ma and 03 denote the gluino mass and 
its phase, and — g^/Air): 

(Amd);,3 = ^e-^('^''+*3)^^Crfm3(m|-mJJ/(m|,m|,m2), 

^ 4v 27rm^^ cos /? ~^ ^' ^ ' j 

(Amrf)^o = -^e-^^'^tan%s<iCrf(m|-m|)^7V*(Ar*tan%-37V;^)m^o 

+ 24.1"jros/3 ^^^^^^^""^^ - 3Ar*)m,o[c^^o(m|,m^o) 

+ slB,{ml,mlo)] (3.30) 



127rm^ cos (3 

In addition, for the 6-quark there is a "pure higgsino" contribution: 

a^mtnib 



^ onm^ sm p cos p 1-1 12 ^ 

X I{ml,ml,mK). (3.31) 



^This convention for the phase of A corresponds to the following normalization of the left-right mixing 
terms in the squark mass matrix: m,|,/j = — m/(>lj + fiRf)- 
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(ii) Corrections to the chcirged lepton masses 



X /(m?^,m?^,mJo) 

+ s2Bo(m2^,mJo)] (3.33) 



4:Trni^ cos /5 



tan 9w X) ^i3^ii'>TT'x° ^e-^o (?7igi , "^^0 ) + CgSo (m~2 , m^o ) 



The two loop functions entering these expressions, / and Bq, are given by 

xy Inx/y + yz\ny/z + xz\nz/x 



I{x,y,z) 



{x - y){y - z){x - z) 



~ , , X In a; — V In v / , v 

B^{x,y)^ 3.34 

x-y 

where the tilde in Bq serves to distinguish this "truncated" version^ of the Passarino- 
Veltman function [39] (following the definitions of [40]). 

Let us briefly discuss these results. Firstly, to gain some intuition on the relative contri- 
butions, it is instructive to consider a limiting case, where the left-right mixing is considered 
a perturbation. The dominant corrections arise from (Am((i,s,b))A3 and (Am6)^+. For the 
first of these, taking ms ~ mj^ ~ mj^, we have (for d — {d, s, b)) 

(^) ^f^!!!|(^tan^5 + A-). (3.35) 

The contribution of these threshold corrections to Sd arises only from Im(/i), as the contri- 
bution from lm{A) is cancelled in constructing — 5^^\ From (3.35), we see that when 
II is negative (^^ = tt) there is destructive interference with the tree-level mass term. 

The presence of an extra diagram sets the b flavor apart from d and s, for which the 
following relation is expected if the soft-breaking sector is flavor-blind: 

d _ — 'nis_^ _ (3.36) 

For the b quark mass, in a regime where the gaugino-higgsino mixing is also small, namely 
/X, rui S> with being the gaugino masses, we also find 

...^ 

^Note that the expression for Bq requires a regulator scale. However, our physical results do not depend 
on this scale. 
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which indicates that when /i is real and positive, so that {Am(^d,s,b))x3 ^Q^^, the contribution 
to 5b is positive when < < tt. 

As is apparent from the hmiting cases above, we observe the appearance of the reparametriza- 
tion invariant CP- violating phases. In particular, at large tan P 



4>a 



(3.38) 



where = Arg (fi), such that the CP-odd phase entering the gluino expression is Arg {firris). 
The analogous behavior of the neutralino contributions is seen most easily in the limit 
/i, rrii » . Then, 

AT* — > 6a e-i^'^""', (3.39) 

and the consequent CP-odd phase in the gaugino piece of the neutralino contribution 
is Arg (/xmi) (or Arg (//y^mirn^))- For the higgsino contribution, the relevant phase is 
Arg (An) = 9a + e^. 

It is worth noting that the origin of the tan /9-enhancement, which is not manifest in 
the general expressions. For the gaugino contributions proportional to the /-function, this 
arises implicitly via the relation 



oc m(itan/3 



(3.40) 



at large tan/?, whereas the mixing angle stays constant in this limit. 

The tan/3 enhancement of the amplitudes involving the Bq function deserves special 
discussion. The relevant diagrams in the mass eigenstate basis appear to be divergent since 
they are proportional to the Passarino-Veltman Bq function. However their tan /3-enhanced 
imaginary parts are finite. As an example, let us consider the chargino contribution which 
involves gaugino-higgsino mixing. The relevant expression appearing in the amplitude is 



1 



cos/3 



(3.41) 



There are in fact several compensating factors of tan (3 in this expression. To sec this, we 
note that Eq.(3.27) implies 



diag(m 



) 



for any (analytic) function /. Therefore, 

1 



cos P 



Bo (rnl,M^+Ml 



21 



(3.42) 



(3.43) 



If Bo here is replaced by a constant independent of M^+, no tan/3 enhancement appears 
since [M^+]2-|^ = V2m^cos(3. This is the reason why the tan /3-enhanced contribution is 
finite and the original Passarino-Veltman function gets replaced by its "truncated" version. 
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It is instructive to expand Bq in powers of {M^+M^+ — rn\^)/m\^ and extract the first 
tan /3-enhanced term. Using the exphcit form of the chargino mass matrix, one finds 

n , /^*^2 ^^tan/5. (3.44) 



This result can be understood via the mass-insertion approximation [41, 42]. To obtain 

tan,/? enhancement, it is necessary to introduce three mass insertions on the chargino line: 
/i - to mix Hi with iJ2) ^^'^ /? - to mix the higgsinos and gauginos, and m2 - to get the 
chirality flip in the diagram. Similar considerations apply to the neutralino contributions. 
We note that analogous calculations have been performed in Refs.[42] and [43]. 

(iii) Scalar-pseudoscaleir mixing 

The computation oiH—A mixing is straightforward, especially in the limit m\ ^ {AH) 
[44] . Non-neghgible contributions arise only from stop, sbottom and stau loops: 



{AH) = 



^ t^ to 



tl t2 
41 ,,|2| /I. |2+„^4 



+ "t'-'lin'^ ^KK) M2e.. + 26,) (3.45) 

m^luPIAT-Ptan^/? 9,9s 
+ -^Z^ m^^2 ^ «/<^ ^i^(2^^, + 29,) 

where the loop function is given by E{a) = — 2 + (a + l)(a — 1)~^ In a. At moderate tan/3 
the stau and sbottom contributions are heavily suppressed due to the factor of in 
the numerator. However, this is not true at large tan /3, as the fourth power of tan $ can 
overcome this suppression. 

It is important to note that {AH) is even under flipping the sign of A, unlike the 
effects of the vertex corrections. Thus, these two sources can interfere constructively or 
destructively, depending on the position in parameter space. For example, choosing 9, = 
and < 9a < leads to destructive interference between Cj;J^^ and C-f^\ 



4 Synopsis of the EDM formulae 

In this section, we compile the different contributions into physical observables measured in 
various EDM experiments. Such observables can be subdivided into three main categories: 
EDMs of paramagnetic atoms and molecules, EDMs of diamagnetic atoms, and the neutron 
EDM. 



4.1 EDMs of paramagnetic atoms — thallium EDM 

Among various paramagnetic systems, the EDM of the thallium atom currently provides 
the best constraints on fundamental CP violation and is often interpreted directly in terms 
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of the EDM of the electron. We therefore speciahze our analysis to the case of Tl, noting 
that the conditions making Cij important for thallium will also be applicable to other 
paramagnetic systems such as the YbF and PbO molecules, where significant experimental 
improvements in EDM measurements arc likely [7, 8]. 

The semileptonic four-fermion operators (2.4) induce the following T-odd nucleon- 
electron interactions^ [9] , 



-eN 



(4.46) 



where the isospin dependence is suppressed. Among these couplings, Cs plays by far the 
most important role for the EDM of paramagnetic atoms because it couples to the spin of 
the electron and is enhanced by the large number of nucleons in heavy atoms. 

A number of atomic calculations [45, 46, 47] (see also Ref. [9] for a more complete list) 
have established the relation between the EDM of thallium, de, and the coefficients of the 
CP-odd electron-nucleon interactions Cs'- 



dTi = -5854 - e 43 GeV x {Cp' - 0.2^^"^). 



(4.47) 



The relevant atomic matrix elements are calculated to within 10 — 20% accuracy [48]. Here 
the relative contribution of the isospin-triplet coupling is suppressed by {N — Z)/Ac:^ 0.2, 
and the effects of Cp and Ct are negligible. 

Before turning to the contributions to Cs, we recall for completeness (and to fix our 
conventions) that the electron EDM de receives a number of well-known contributions from 
superpartner loops. At one-loop order, threshold corrections from e — x"^ and e — x° loops 
dominate, and we follow the notation of Ibrahim and Nath [14] for these terms (up to a 
different convention for the squark mixing angles as specified earlier). At two-loop order, 
non-negligible Barr-Zce-type Higgs-mediated graphs contribute [49], where CP violation 
enters through the couplings of (s)fermions to the pseudoscalar Higgs A. We follow the 
work of Chang, Keung and Pilaftsis [49] (though our phase conventions are somewhat 
different) and to emphasize an additional contribution, we quote the full expression below: 



a 



rrie 



F 



31 



m-. 



32 



tan /3, 



(4.48) 



where j — t,b,T are flavors running in the sfermion loop. Nf — Ni, — 3, and Nt- — 1, 
accounts for the trace over color, while the two-loop function F{z) is given explicitly in 
[49]. The CP- violating couplings in our conventions are: 







sm(2ej)mjRjIm{iie-''l'i) 
sin /3 cos /3 



(4.49) 



Note that while there are in principle a large number of possible diagrams in this class 
[50] , those containing stop, sbottom, and stau loops should be dominant at large tan /3. We 

^Note that we differ from Ref. [9] in the definition of 75, so that we use Pj, = (1 — 75)72. 
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emphasize here that the stau-loop contribution, which has hitherto been neglected in the 
hterature, is a significant addition - in fact for the CMSSM and NUHM at tan /3 = 50 to 
be considered here, it often provides the dominant contribution due to the generically hght 
staus. 

In order to see the dependence on the reparametrization invariant phases, it is useful to 
again consider the small- mixing regime where m| ~ rrij^ ^ m\, m\ji. For large tan /3, the 
EDM then reduces to 



d'r-'^^^ij^tM - -Q^mrf'^^^^^^^ I 7:1 1 sin(e^, +e,)tan"/3, (4.50) 




m 



where wc Tised — —1, and defined n as n = 1 for j — t and n = 3 for j = b,T. It is 
important to recognize that when = the two-loop contribution to the electron EDM is 
negative which interferes destructively with the one-loop contribution referred to above^ . 

Having accounted for de, the calculation of Cs induced by (2.4) follows from the low- 
energy theorem that expresses the nucleon matrix element of agG^iyC^" in terms of the 
mass of the nucleon and the QCD beta function [52] . As a consequence of this theorem, for 
a heavy quark Q, {N\mQQQ\N) ~ 70 MeV. This theorem is valid in the exact chiral limit, 
and thus receives corrections from the non-zero values of the light quark masses. Following 
[26], we parametrize the m^ss matrix element as k = (A^|msSs|A^)/220 MeV. Using this 
parametrization, as well as the measured value of (m„ + md){N\uu + dd\N)/2 ~ 45 MeV 
and the ratio of the quark masses rriu/md ~ 0.55, we arrive at the following general formula 
for Cs in terms of Cge, 

_ 29MeV K X 220 MeV , ^ 66 MeV(l - 0.25k) 
rud rus mi, 

ritrip ^ 1-3 MeV 

Cs ~ -Cde- 



rrid 



In this formula, Cdejvad and Cse/Tris should be taken at the normalization scale 1 GeV, while 
Cbe/mb must be normalized at nib- In the case of minimal SUSY models, these combinations 
are scale invariant, as Cqe ~ rriq. There is a rather large uncertainty in the isospin triplet 
coupling due to the poorly known value of (m^ — mu){p\uu — dd\p). Fortunately, C^"^ 
is numerically small and a further suppression via (4.47) makes the contribution of Cg^*^ 
completely negligible for all applications. Thus, from now on, we neglect C5 and drop 
the superscript in C^"^. 

The largest uncertainty in C5 (4.51) that may significantly affect numerical values of 
EDMs originates from the poorly known value of k. The assumption that the s quark 
behaves as a heavy quark would give k, ~ 0.3, which we regard as a lower bound. The 
analysis of the baryon octet mass splitting at leading order in chiral perturbation theory 
instead suggests k ~ 1, while an improved next-to- leading order calculation gives a predic- 
tion K, ~ 0.50 ± 0.25 [53] that we use as the central point for our numerical analysis. From 

^Note that Barr-Zee-type diagrams with chargino loops also contribute for ^ [51, 37]. 
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these numbers it is clear that the overall precision with which this matrix element can be 
estimated is not better than 30%. 

Substituting Eqs. (3.22) and (3.24) into (4.51), we arrive at the final expression for Cs 
in terms of m^, Am^, 5i and {AH): 



5.5 X l0-i"+->^^2 



taii2/3 



m^|l + Anie/me 



(O)l 



(1 - 0.25/x) 
|1 + Amfe/m^' 



(0)| 



+ 



|1 + Anis/ms 
0.5 



(o)i 



sin(5s - 6e) 



{AH)\ 



{AH)\ 



(4.52) 



|1 + Amd/m\ 



(O)l 



m\ ) 



We use this calculation of C5, along with the standard SUSY calculations of reviewed 
above, in the analysis of the thallium EDM as a function of SUSY masses and phases which 
will be discussed in Section 5. 



4.2 Neutron EDM 

The calculation of the neutron EDM in terms of the Wilson coefficients of Eq. (2.5) repre- 
sents a difficult non-perturbative problem. The most common approach invokes "naive 
dimensional analysis" (NDA) [54, 55] which gives an order of magnitude estimate for 
dn{di, di, w) with uncertain signs. Some partial calculations have also been performed with 
the use of chiral perturbation theory [56, 57] and lattice QCD [58]. Recently, a universal 
treatment of all operators within the same method was developed in Refs. [59, 60] using 
QCD sum rule techniques. We briefly summarize these results below. 

It is natural to expect that the dominant contribution to the neutron EDM comes from 
the EDMs and color EDMs of u and d quarks. We recall that the ^-dependence is removed 
by the Peccei-Quinn mechanism, which at the same time reduces the contribution of the 
color EDM of the s-quark [59]. The QCD sum rule analysis [59] leads to the following 
result: 

dn{di, dg) = 0.7{dd - 0.25rf„) + 0.55e(Jrf + 0.5J„), (4.53) 

where the value of the quark vacuum condensate (qq) = (225 MeV)^ has been used. Here 
dq and dq are to be normalized at the hadronic scale which we assume to be 1 GeV. The 
relation to the Wilson coefficients at the weak scale is as follows, 

dq{l GeV) ~ 0.91dq{Mz); dq{l GeV) ~ 1.2dq{Mz). (4.54) 

Note that dq as defined by Eq. (2.2) has a much milder QCD scaling than an alternative 
definition for dq which includes a factor of Qs in its definition. The reader should also note 
that the quark masses used for the SUSY calculations of dq and dq should be taken at the 
weak scale, where their numerical values are smaller than the low energy values by a factor 
of ~ 0.35, e.g. md{Mz) ~ 9.5 MeV x 0.35. 
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The expression (4.53) has several nice features. Its flavor composition agrees with the 
predictions of the SU(6) constituent quark model, while the proportionality to dq{qq) ~ 
''^qifil) ~ /^"^^ removes any sensitivity to the absolute value of the light quark mass. The 
quark EDM part of (4.53) is in good agreement with lattice calculations [58]. Finally, the 
overall uncertainty of (4.53) is estimated to be at the 30%-50% level [59]. For completeness, 
we recall that there are standard SUSY one- and two-loop contributions to dq and dq in 
analogy with those discussed for the electron EDM above. For 1-loop quark EDMs and 
CEDMs we follow [14] and note that, in comparison to the electron EDM, for quarks there 
are additional one-loop diagrams containing squark-gluino loops. The two-loop Barr-Zee- 
type contributions to the quark EDMs [49] also receive a significant correction from the 
stau-loop in addition to those containing stops and sbottoms. 

The Weinberg operator can also be an important source of CP violation, especially when 
the third generation of squarks is lighter than the first two generations. Unfortunately, a 
reliable calculation of its contribution to the neutron EDM is problematic at this point. 
The QCD sum rule approach [60] produces an estimate that is close to the NDA result but 
this calculation is at the border-line of applicability of the method. With w normalized at 
the 1 GeV scale one finds: 

dn{w) ~ 20MeV xew. (4.55) 

This estimate is assessed to be vahd within a factor of 2-3 [60] . However, it seems unhkely 
that this uncertainty can be significantly reduced. An important implication is then that 
even the sign cannot be reliably inferred, as the next order corrections to the QCD sum 
rule determination (4.55) are not calculable and plausibly are large. 

When tan j3 is large, we expect that a generically dominant contribution to the Weinberg 
operator will be induced by the 6-quark color EDM (see also [61, 62]). Thus, the expression 
for w{db) linearly enhanced by tan (5 takes the form: 

w{l GeV) = 0.72«;(m,) = -0.72 x = -0.68 x (4.56) 

327r^mf, 327r^m{, 

It is well known that the t and c quark two-loop contributions are just a fraction of the total 
EDM at low tan/3 [63]. Since they do not grow with tan/3, their effect is expected to be 
negligible but nonetheless they are included in our analysis. We remark at this point that 
it is common in the SUSY-EDM literature to quote the Weinberg operator at the hadronic 
scale as 3.3 x w{Mz)- This relation originates from Ref. [62] where the contributions to w 
induced by and dc were added at the h and c thresholds with some specific assumptions 
about the SUSY spectrum and a fixed (rather low) value of tan /3. For any other point 
in the parameter space this coefficient is going to be different, and therefore an explicit 
summation of the h and c thresholds should be performed separately. 

There are many different types of four-fermion operators generated via the vertex cor- 
rection mechanism that contribute to the neutron EDM at leading order: Cds^ Csd, Cm, Cat, 
Cbs, and Csb- The contribution of the first two operators to dn could be important a priori 
[35]. We note, however, that the CP- violating phases induced under renormalization are 
the same for the d and s quarks to very high accuracy within the framework of a fiavor-blind 
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soft-breaking sector. Thus, Ss ~ Sd and Cds ~ —Cgd ~ 0, from Eq. (3.36). This however 
does not apply to Cm and Cdb or C^s and Csb due to an additional chargino contribution 
to the 6-quark mass. The estimates of dn induced by Cm and Cat can be obtained along 
the lines of Ref. [60]. We supply some of the details of this calculation in Appendix B, and 
quote here the result: 

, , e 0.65 X lO-'^GcV^^ 

4 Cm Cm, 4.57 

nib 

where we took into account that Cm — ^Cdb- The Yukawa couplings Yd and Yb that enter 
the coefficient Cbd in (4.57) are normalized at 1 GeV and uib respectively. In the final result 
for dni Cbs and Csb may be as important as Cbd and Cdb- However, we are not aware of any 
reliable way to estimate dn{Cbs, Csb)- 

Combining the different contributions, (4.53), (4.55) and (4.57), we arrive at a final 
formula for the neutron EDM as a function of the different Wilson coefficients: 

dn = dn{dq, dg) + dn{w) + dn{Cbd)- (4.58) 

One should keep in mind at this point that, as emphasized above, Eqs. (4.55) and (4.57) are 
undoubtedly of poorer precision than (4.53). Therefore, a reliable calculation of dn in terms 
of a specific combination of SUSY CP-violating phases is possible only if both the second 
and third terms in (4.58) are smaller than dn{dq,dg). When instead they contribute at a 
level comparable to the EDMs and color EDMs of quarks, an interpretation of the neutron 
EDM bound in terms of constraints on specific CP-violating phases becomes problematic. 



4.3 EDMs of diamagnetic atoms — mercury EDM 

EDMs of diamagnetic atoms, i.e. atoms with total electron angular momentum equal to 
zero, also provide an important test of CP violation [9]. The current limit on the EDM of 
mercury [3] furnishes one of the most sensitive constraints on SUSY CP-violating phases 
[16]. However, the calculation of dug is undoubtedly the most difficult as it requires QCD, 
nuclear, and also atomic input. 

The atomic EDM of mercury arises from several important sources, namely, the Schiff 
moment S [64], the electron EDM d^. and also the electron- nucleus interactions Cs and Cp 
(see, e.g. Ref. [9] for a comprehensive review). Schematically, the mercury EDM can be 
represented as 

dug = dug (^S[gT:NN {di, C,,,J], Cs[Cqei Cp[Ce,], 4) , (4.59) 

where g-jrNN collectively denotes the CP-odd pion-nucleon couplings. The atomic and nu- 
clear parts of the calculation have been performed by different groups, and several results 
such as dugiS) [65] and S{g.j^NN) [66] have been updated recently. The most important 
numerical change comes from a new QCD sum-rule calculation of g-wNN{di) [67], that ob- 
tained a preferred range and "best" value for this coupling. Previous estimates [33, 16] are 
within this preferred range, but close to the largest possible value for the g-KNN{di) matrix 
element. The best value that follows from the sum- rule analysis [67] is a factor of ~ 2.5 — 3 
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smaller than previously used values [16]. Incorporating the changes in the atomic matrix 
elements, the current theoretical estimate for the mercury EDM induced by dimension 5 
operators stands as: 

c^Hg = 7 X 10"^e(J„ - dd) + 10"^ X 4. (4.60) 

Note that the electron EDM also contributes to dug, although less significantly than to 
d-n which provides a more stringent bound. Therefore, the most valuable feature of dng is 
its sensitivity to the triplet combination of color EDM operators 4, which surpasses the 
neutron EDM sensitivity to this combination of these operators by a factor of a few. The 
overall uncertainty in the QCD part of the calculation is considerably larger than that for 
the neutron due to significant cancellations, refiecting the fact that the result vanishes in the 
vacuum factorization approximation (see Appendix B). However, the dominant dependence 
on the {du — dd) combination ensures that these uncertainties enter as an overall factor and 
therefore do not significantly alter the shape of the unconstrained band of the parameter 
space in the 9^ — 9a plane. 

One should note that the mercury EDM also receives contributions from the color EDM 
of the strange quark [16] and the EDMs of light quarks, although their contribution is 
subleading. Finally, the Weinberg operator does not provide any appreciable contribution 
to (ing because its contribution to Qt^nn is suppressed by an additional chiral factor of 
mjl GeV ~ 10-2. 

At large tan/3 one should also include additional contributions to the Schiff moment 
coming from the four-quark operators, as well as the effects of the semi-leptonic operators 
Cs and Cp: 

dug = 7 X 10-2 e {du - dd) + 10-2 de 

-1.4 X lO-^eGeV^ + 3.3k^ + ^(1 - 0.25k)) (4.61) 

V rrid rris rrih J 

+3.5 X lO-^QeYeCs + 4 x lO^^GeV eCp. 

The second line in (4.61) is the contribution of the four-quark operators to the Schiff 
moment. The details of the g-nNN^Cq^q^) estimates, based on a factorization assumption, are 
given in Appendix B. The result depends on the same parameter k (the nucleon strangeness 
in the 0+ channel) as the coefficient Cg- As emphasized before, Cgd and Cdd are not induced 
via vertex corrections, and arise only through A — H mixing. Finally, the contribution of 
Cp is given primarily by C^d [26] while the contributions from C^s and C^b are consistent 
with zero. The QCD normalization of all operators is the same as in the thallium and 
neutron cases. 

In the region of parameter space where the color EDM contribution to dug is dominant, 
the mercury EDM hmit is extremely valuable since, as noted above, it constrains the triplet 
combination of up and down color EDMs, and therefore a calculable combination of the CP- 
violating phases. On the other hand, when other terms in Eq. (4.61) contribute at a level 
comparable to that of 4 — 4, which in fact is expected at large tan/3, different uncertainties 
no longer factorize in front of a given phase combination, and the interpretation of the 
mercury EDM constraint becomes more difficult. 
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5 Numerical analysis of EDMs at large tan j3 



In this section we analyze all three observables, dxi, dn and rfng, within two classes of 
MSSM models. In both cases we assume flavor-blind SUSY breaking, a common trilinear 
soft-breaking parameter Aq, and a universal (real) gaugino mass parameter mi/2 at the 
GUT scale. This narrows down the number of SUSY CP- violating phases to two, which 
can be identified with the phases of the /i and Aq parameters. We then perform an EDM 
analysis separately for these two phases. In all numerical runs the value of tan/3 was chosen 
to be nearly maximal, tan (3 = 50. We defer an analysis of the constraints on the SUSY CP- 
violating phases to a subsequent publication [68] and now focus on the relative importance 
of various EDM contributions. 

The constrained MSSM (CMSSM) is a popular framework defined by the following set 
of universal SUSY parameters at the GUT scale: 

{tan/3, mo, mi/2, \Ao\, Oa, 6'^}- (5.62) 

The magnitude of the /x-paramcter and the pseudoscalar mass are determined by the ra- 
diative electroweak symmetry breaking conditions (i.e. by reproducing the observed value 
of Mz). We note that, in models with CP violation, the usual sign ambiguity of becomes 
a phase ambiguity. 

The renormalization group running from the GUT scale introduces considerable mass 
splittings in the spectrum of superpartners. For example, the gluino becomes much heavier 
than (roughly triple) the rest of the gauginos. The RG running can also make the scalar 
quarks quite heavy, 

ml^{Mz) ~ + 6ml/, + 0{Ml), (5.63) 

especially if mi/2 is large. In our analysis, the masses of the A and H Higgs bosons are 
functions of the above input parameters and are of particular importance. In the CMSSM 
with large tan p, they tend to be as heavy as the sleptons. As a consequence, the parameter 
^ determining the relative importance of the four-fermion operators (2.6) cannot be varied 
over a large range once tan /3 is fixed. 

In contrast, ^ becomes a free parameter in another version of the MSSM, a CMSSM-type 
model with non-universal Higgs masses (NUHM). In this model the squarks and sleptons 
are still mass degenerate at the CUT scale, but the Higgs soft masses are not. The model 
differs from the CMSSM in having two additional Higgs mass parameters, which can be 
taken to be the mass of the pseudoscalar and /i evaluated at the weak scale, so that 
the parameter set becomes 

{tan/3, mo, my,, tjia, |^o|, Oa, 6*^} ■ (5.64) 

Relaxing the scalar mass universality brings an interesting new degree of freedom into 
weak-scale phenomenology while preserving the flavor-degeneracy of the CMSSM. For this 
study it is important that msiepton/nT'A and msquark/^A can be varied over a large range. 
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thus presenting an opportunity to study the transition to a regime of dominance for the 
four-fermion operators. 

In what follows, we present our numerical results for the EDMs using three two- 
dimensional slices through the parameter space - planes of mi/2-mo (typical of the CMSSM) 
and /x-mA and mA-niQ (as examples of the NUHM). For each EDM observable, we plot one 
CMSSM figure with = and Oa = 7r/2. Note that the phase of A runs in the RGEs, so 
that at low energies, the phases of the various Ai are different and those of the squarks are 
typically <^ 7r/2, whereas the phases of Aj remain relatively large. In contrast, the phase 
of n does not run so that 9f^{EW) = d^{GUT). We also plot for each EDM, three different 
NUHM cases: one for 6*^ = 10~^7r and 6*^ = and two different parameter planes for 6'^ = 
and 9 A = 7r/2. In all cases we fix tan /? = 50 and \Aq\ = 300 GeV as well as rrit = 175 GeV 
for the pole mass, and = 4.25 GeV for the running MS mass evaluated at itself. 
Note that our results for the EDMs are not particularly sensitive to the choice of \Aq\. 

In all cases, we apply a set of phenomenological constraints that we describe briefly. 
We apply the constraints on new particles from direct LEP searches, namely m^± > 
104 GeV [69] and Higgs mass limits [10] of irih > 114 GeV. The light Higgs mass has 
been computed using FeynHiggs [70]. We also require the calculated branching ratio for 
6 — > S7 to be consistent with experimental measurements [71]. For each set of parameters, 
we calculate the relic density of neutralinos having set the CP-violating phases to 0. As 
such the domains we display for a given relic density are only indicitivc. However, we note 
that in the CMSSM (and more generally when sfcrmion masses are not too degenerate 
at the weak scale), the relic density of binos is not significantly affected by these phases 
[13]. For the relic density of neutralinos x, we show regions of the parameter planes with 
0.1 < fl^h"^ < 0.3. We use this more conservative range in addition to that suggested by 
the recent WMAP data [72] (0.094 < n^h'^ < 0.129) for clarity in the figures. The effect of 
the WMAP densities in the CMSSM was discussed in [73] and in the NUHM in [74]. We 
also require that the lightest supersymmetric particle (LSP) be a neutralino. For further 
details on these constraints see [75] . 

5.1 Thallium EDM 

As emphasized in section 4.1, the thallium EDM receives two main contributions at large 
tan/?, from dg and Cs, while the EDM of the electron itself can be generated in different 
ways at both one- and two-loop order. Figure 4a plots contours of constant dT\{Cs)/dTi 
(sohd, blue) and (dashed, red) in the CMSSM as a function of mi/2 and mo. 

The GUT-scale phase of is set to be maximal, 9a — 7r/2, while 9^ — 0. The dark (red) 
shaded lower-right half of the plane corresponds to a stau-LSP region which is excluded if 
the LSP is stable. The medium (green) shaded region is excluded by 6 57 (computed 
with non-zero phases). In the broad light shaded region the relic density of neutralinos is 
0.1 < Q,h^ < 0.3 (computed in the absence of phases) and the thin slightly darker region 
shows the WMAP relic density. In this figure, the rehc density is determined primarily by 
co-annihilations with the nearly degenerate f and/or strong s-channel annihilation through 
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Figure 4: Contributions of four-fermion operators (C5) and the electron EDM (de) to the 
EDM of thalhum as functions of the SUSY parameters. Blue (solid) lines are contours of 
constant d'T\{Cs)/dTi, and red (dashed) lines are the contours of constant d1~^°°^ / d^°^"-K 
A thick solid black line marks the current experimental constraint on \di\\, and black 
dashed and dot-dashed lines represent constraints from direct searches for charginos and 
the Higgs. The legend indicates the color(grey-scale)-coding of the excluded regions and 
preferred domains as described in the text. 
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the heavy Higgses H, A. In almost all of the unshaded portions of this figure, the rehc 
density is too high and these parameter choices would be excluded if the LSP is stable. 
Also shown by the nearly vertical black dashed line is the LEP bound on the chargino mass 
(smaller values of mi/2 are excluded) and the black dot-dashed curve is the LEP bound on 
the Higgs mass (which also exlcudes small values of mi/2 and mo). 

The importance of Cs and dl~''°°P grow as one moves to lower mi/2 and higher mo as 
indicated by the blue (solid) and red (dashed) lines. In particular, the contribution from 
Cs changes dramatically, from merely a 5% effect to a 100% contribution, while the two- 
loop contributions are generically large. Moreover, since the effects of Cs and d1~^°°^ have 
the opposite sign relative to the one-loop contribution to dg, the overall consequence is a 
significant reduction of the predicted observable thallium EDM over a large portion of the 
parameter space. In fact, the total EDM of thalhum is generically several times smaller 
than either di\{Cs) or d^\{d^°*°'^) individually. The band in parameter space where there 
is greater than 80% cancellation between Cs and d^ is shaded in blue (dark grey), i.e., 
dT:\{C s) / dfx > 5. In the center of this band, rfxi actually vanishes. Our results for Cs are 
weakly dependent on k. The 50% uncertainty in k, translates into a 10% uncertainty in 
Cs- This lack of sensitivity to k, can be inferred from Eq. (4.52) where the first term in the 
square bracket is dominant. 

As noted earlier, the importance of the two-loop d^ contribution stems primarily from 
stau loops since the staus are generically lighter than the rest of the squarks and sleptons. 
The left-most band, shaded green (light grey), in the parameter space shows the region 
where these two-loop contributions nearly screen the one-loop effects in de (to 95%). In 
this band, dl'^""^ /de > 20. In the center of the band, de — ^ and the observable atomic 
EDM is given entirely by Cs- Note that Cs is dominated by vertex corrections, C\Y\ and 
the contribution due to H-A mixing is not very important. 

Turning now to the NUHM, we observe that the cosmologically preferred region of 
parameter space can be substantially wider. Figure 4b presents the thallium EDM results 
in the /i—mA plane, and one can see that even for uia as high as 500 GeV Cs may contribute 
10% of the total thalhum EDM, while for uia < 300 GeV Cs often provides the dominant 
contribution. Here, we have fixed mi/2 = mo = 500 GeV. For this choice of parameters, 
positivity of the Higgs soft masses at the GUT scale would require /i < 1600 GeV [30]. In 
this figure, the r-LSP region is now in the upper left corner and the region excluded by 
6 — > S7 extends down from that. The region with the relic density between 0.1 and 0.3 
are the two funnel-like strips which become nearly horizontal at large /i- This region is 
allowed by the rapid s-channel annihilation of x's through heavy H, A. Note that ^ 
0.4mi/2 ~ here. Between these two strips, the relic density is less than 0.1 and so 

is not excluded. The WMAP region corresponds to the thin strips in the interior of the 
funnel. Once again, we see the blue band with ciTi(C's)/c^Ti > 5 at uia ^ 250 GeV which 
now intersects the WMAP region at /i ~ 400 GeV and m^ ~ 200 GeV. At slightly lower 
niA, the 2-loop contribution cancels the 1-loop contribution to de and dxi = dT\{Cs)- 

Perhaps the most interesting slice to consider is the ttia — mo plane plotted in Fig. 4c. 
From Fig. 4b, we see that in order to maximize the size of the WMAP cosmological region. 
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we should choose = 400 GeV. Here, mi/2 is fixed to 600 GeV. Indeed, in Figs 4c and 
4d, we sec a large region with ttia ^ 300 GeV (and ~ 650 GeV) where flx^"^ ~ 0.11. 
For intermediate values of the relic density is small and a dark matter candidate is 
not provided. The f-LSP region now excludes the lower portion of the figure and 6 — > 57 
excludes a small area at higher mo when is large. 

In Fig. 4c, one can see that over most of the cosmologically preferred region, Cs provides 
the dominant contribution to the atomic EDM. It follows that nowhere within this region 

can the observable EDM be interpreted purely in terms of d^. alone. Moreover, one observes 
that within the same cosmologically preferred region there is a large band of the parameter 
space where the two contributions to the total atomic EDM cancel, implying no constraint 
on the CP-violating phase. Once again, 50% changes in k lead to 10% changes in Cs- 
The final figure (4d) presents the same shce of parameter space within the NUHM as 
Fig. 4c for = and = 10~^7r. In this regime, the one-loop contribution to de is 
tan /3-enhanced, and provides the dominant effect. We see that even for ~ O{10^^) the 
domain of parameter space with mo < 1100 GeV is excluded by the thallium EDM bound 
as shown by the thick solid black curve. The contribution of Cs reaches 20% on this plot 
for low rriA and high values of mo. Glearly, in this case, over the range of parameter space 
that is potentially accessible at future colliders, the one-loop electron EDM provides the 
leading contribution to the atomic EDM of thallium. Here, we note that when 7^ 0, the 
sensitivity to k, is larger and a 50% change in k, leads to a comparable change in Cs- 

In summary, this numerical analysis indicates that although Cs plays a very important 
role for paramagnetic EDMs at large tan /3, the precise values are somewhat (a factor of a 
few) smaller than one would naively expect from the simple scaling arguments presented 
in Sec. 2. The main reason for this is that Cs is driven mostly by squark loops while de 
also receives contributions from selectrons, and as a consequence Cs tends to be further 
suppressed due to the generically heavy masses of the squarks. Another important factor 
is the renormalization group evolution from the GUT scale down to the electroweak scale 
that considerably suppresses the size of Im^l^ relative to ImAg leading to an additional 
suppression of Cs/d^. 

5.2 Neutron EDM 

Repeating the same sequence of plots for the neutron EDM, we observe important differ- 
ences with respect to thallium. In all four panels of Fig. 5, the phenomenological and 
cosmological considerations are identical to those described for Fig. 4. In the GMSSM, 
with 9a — 7r/2, the current experimental bound, shown as the thick black curve, for dn 
excludes much more of the plane than that for dxi (due to the absence of any significant 
cancellations). Indeed, for low values of mi/2 around 100 GeV, scales for mo above 1 TeV 
are probed by the experimental constraint. In the other planes, the choices for mi/2 and 
mo and/or were made so as to satisfy ab initio the current experimenatal bounds over the 
displayed plane. Furthermore, there is httle sensitivity to 9^ — lO'^n compared to that of 
the thallium EDM. The contribution of the two-loop quark EDMs and GEDMs is in gen- 
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Figure 5: Contributions of the Weinberg operator, dn{w) = and four-fermion opera- 
tors, dn{Cij) denoted as Cbd here, to the neutron EDM as functions of SUSY parameters. 
Blue (solid grey) lines denote contours of constant dn{d'^~^"°^)/dn, while red (dashed) lines 
are contours of constant dn{Cbd)/dn, and finally green (dot-dashed) lines are contours of 
constant dyj/dn- The thick solid black line delineates the experimental bound on and 
black dashed and dot-dashed lines represent constraints from direct searches for neutrali- 
nos and the Higgs. 
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eral very small, not exceeding 10% of the total Moreover, of particular relevance here, 
a scan over the parameter space of the CMSSM and NUHM reveals that the four-quark 
interactions and the Weinberg operator contribute at the 10% level at most. Both the 
Weinberg operator and the two loop contributions df~''°°^ scale as (two-loop factor) x tan /3 
and therefore it is natural that they both contribute at a similar level. The contribution of 
four-fermion operators is dominated by C^d and increases when mo is increased and/or 
is lowered. Although outside the parameter regime included in the plots, it is worth noting 
that a factor of 3 increase in the maximum value of mo to around 3 TeV in Fig. 5c (at low 
ttia) would bring the Cm contribution to a level comparable with dn{di,di). 

We conclude that for all mo < 1 TeV the neutron EDM at large tan f3 can be reliably 
calculated in terms of the the EDMs and CEDMS of quarks. The hadronic uncertainties 
associated with the matrix elements of four-fermion and Weinberg operators, which we can 
only estimate within a factor of 3, do not significantly disrupt the predictions for dn- This 
means that even at large tan f3 the calculation of dn can provide meaningful constraints on 
specific combinations of the quark EDMs and CEDMs. 

5.3 Mercury EDM 

In many ways, the mercury EDM occupies an intermediate position between dn and dxi- 
Results for the Hg EDM are shown in Fig. 6. Once again, the phenomenological and 
cosmological considerations are identical to those described for Fig. 4. In the CMSSM with 
6a ^ and = 0, the current experimental constraint on d^g shows even greater sensitivity 
to mo than does dn- On the other hand, the contribution of four-fermion operators to dug 
is relatively less important than for dxi, but more important than for d„. It is interesting 
to note that the contribution of Cs has a similar magnitude, but the opposite sign, to that 
of C<jiy2) and over a large part of the parameter space these two effects tend to compensate 
each other. An exact cancellation would be an accident, as Cg^g^ has considerable hadronic 
uncertainties associated with the calculation of the pion-nucleon CP- violating coupling and 
the Schiff moment. The contribution from Cp is negligible. 

The effects of four-fermion operators become much more pronounced in the NUHM 
model. Fig. 6c shows that for low m^ the mercury EDM receives contributions of up to 50% 
from Cs and Cq-^q^ (for the chosen parameters (iHg(C's) ~ ~'^Hg(C'(?i(?2) and only Cs is shown 
on the plot). These contributions drop quickly with increasing rriA and for m^ > 250 GeV 
do not exceed 10% of the total EDM. It is also important to note that for both the CMSSM 
and the NUHM the electron EDM remains subdominant for 9a ^ and = 0. This 
provides us with some confidence in the calculation of the mercury EDM at 9a 0, ^/i = 
in the range niA > 250 GeV and mo < 1.5 TeV where both d^. and four-fermion operators 
provide contributions of less than 10%. We conclude that in this part of the parameter 
space the mercury EDM is dominated by contributions from [d^ — dd)- 

This situation changes when we consider the case of 9a = 0, 9^ ^ 0, shown in Fig. 
6d. One can readily see that the electron EDM becomes important everywhere within the 
domain of parameter space included in the plot, and its inclusion leads to a significant 
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Figure 6: Contributions of the four-fermion operators and the electron EDM to the 
EDM of mercury as functions of the SUSY parameters. Blue (solid grey) lines denote 
contours of constant d}ig{Cs)/du_g, while red (dashed) lines are the contours of constant 
'^Hg('^e)/'^Hg (indicated as de/dug), and green (dot-dashed) lines are contours of constant 
dHgiCqiq2) / dug (indicated as Cbd/dug on the plot). The thick solid black line denotes the 
experimental bound on {dugl, while black dashed and dot-dashed lines represent constraints 
from direct searches for neutralinos and the Higgs. 
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reduction in the total mercury EDM. In practice, however, in this regime c^ti is dominated 
by df, and, on taking the bound on d^i into account, the effect of rfg on rfng will be suppressed. 
Wc note finally that, unlike the previous observables, at 0^ = 0, 6^ ^ the mercury EDM 
also depends sensitively on Cs and the four-fermion contributions to the Schiff moment. 

6 Conclusions 

In this paper, we have presented the results of SUSY-EDM calculations in the large tan (3 
regime. Large values of tan/3 are motivated by a variety of theoretical models and the 
SUSY Higgs searches at LEP also suggest that tan/3 ^ 5. Our main goal was to combine 
the existing one and two-loop calculations of the EDMs of the constituents (electrons, 
quarks, etc.) with the effects of four-fermion CP-odd interactions, the importance of which 
was recently emphasized in Ref. [26]. We concentrated on comparing the relative size of 
these contributions, while a more detailed analysis of the constraints on SUSY phases will 
appear elsewhere [68]. 

We have compiled the most accurate formulae to date to obtain the three observables: 
the EDMs of the neutron, mercury and thallium as functions of the Wilson coefficients, 
which in turn are expressed in terms of SUSY parameters. Mostly these are compilations 
of existing QCD/nuclear/atomic calculations, although in certain cases such as dn{Cbd) and 
dng{Cbd), new QCD estimates have been obtained. On the SUSY side, we have performed 
a complete calculation of 1-loop corrections to the quark and lepton masses in the large 
tan P regime, including the effects of SUSY CP- violating phases. These mass corrections 
(or rather their imaginary parts) determine the strength of the CP-violating Higgs-fermion 
interaction, the leading source for the CP-odd quark-quark and quark-electron interactions. 
At the subleading level, only the effects of scalar-pseudoscalar mixing were taken into 
account as these are important for low pseudoscalar masses niA, although a large class 
of box diagrams can contribute as well. We also included the stau-loop contribution to 
the two-loop corrections to the EDMs, and recalculated the corresponding sign to firmly 
establish the pattern of interference between one- and two-loops. 

The impact of the four-fermion operators is most striking for the thallium EDM. We 
observe that for low and high mo the electron EDM is no longer dominant, and Cs 
provides as large a contribution (or possibly even larger) to d^i as does dg. A similar conclu- 
sion will apply to other paramagnetic species as well. This is an important point because it 
illustrates that in SUSY models with large tan f3, the EDM of a heavy paramagnetic atom 
or molecule cannot be entirely attributed to d^.. Should the next generation of paramagnetic 
EDM experiments [7, 8] detect a non-zero signal, only a combination of positive measure- 
ments with different atoms or molecules could unambiguously separate the two effects. It 
is also noteworthy that the total electron EDM is reduced due to destructive interference 
between the one- and two-loop contributions, particularly due to f-loops. 

In the Cg-dominated regime for dxi, the main uncertainty in the theoretical calculation 
comes from the uncertainty of the strange quark scalar matrix element over the nucleon. 
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This translates into a 10% uncertainty when the EDM is induced by 9a and 30 — 50% 
uncertainty for the ^^-induced EDM. The transition to dominance of for Q occurs at 
a very large mass scale for the superpartners, and in reality for uiq < 2 TeV the contribution 
of Cs{9fj) is always subdominant. Yet, in the case that the EDM is induced by 9a, we find 
that the four-fermion operators can easily contribute as much as the electron EDM^. We also 
find that, in certain parts of the parameter space, the physical observable dxi is considerably 
smaller than the individual contributions. 

In the case of the mercury EDM, there are a variety of different contributions that have 
to be taken into account. We find that at large tan (3 the contribution of Cs and other four- 
fermion operators may become important. From our analysis it is clear that the main prob- 
lem with the interpretation of the mercury EDM at large tan /? is that a rather large number 
of contributions appear with many of them having different hadronic/nuclear/atomic ma- 
trix elements. Therefore, a reliable calculation (i.e. better than 100% uncertainty) becomes 
problematic at large tan/3. This contrasts with low and moderate tan/9 where the EDM 
of mercury is usually dominated by the color EDMs of the light quarks, and the hadronic 
uncertainties are a less harmful overall factor. Nevertheless, for 9^ = 0, uia ^ 250 GeV and 
mo ^ 1.5 TeV the effects of four-fermion operators and the electron EDM are subdominant 
and a meaningful estimate for dug is possible. 

We determined that the EDM of the neutron is the least susceptible observable, with 
the contributions of four fermion operators never exceeding 10% of the total EDM for 
niQ < 1 TeV. On the other hand, for 9a ^ the contribution of the Weinberg operator 
can be quite substantial. According to our estimates for the matrix elements, it is < 
10% of the total EDM in much of the parameter space considered (if an NDA estimate 
of dn{w) is employed, this contribution is larger). Although subleading, we infer that at 
large tan/3 the calculation of the neutron EDM can be prone to uncertainties related to 
the matrix element of the Weinberg operator. It is also worth noting that the dominant, 
tan /3-enhanced, contribution to the Weinberg operator comes from the color EDM of the 
6-quark. 

Finally, we remark that low values of ^ 300 GeV can be preferred by considerations 
of the cosmological neutralino abundance. This is related to the new annihilation channels 

due to Higgs exchange that open up when niA is not excessively large. In that part of 
the parameter space the effects of four-fermion operators are very important, and therefore 
models that include both CP violation and neutralino dark matter should incorporate the 
effects of four-fermion interactions into the EDM constraints. 
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Appendices 

A Neutron EDM induced by four-fermion operators 

As discussed in Section 4, the direct QCD sum-rule analysis for the neutron EDM becomes 
problematic for sources of dimension six, due in part to the presence of unknown conden- 
sates. For the contribution of the Weinberg operator, one can nonetheless obtain a useful 
estimate [60] by relating dn to the neutron anomalous magnetic moment, /in, by a chiral 
rotation of the nucleon wavefunction. For similar reasons, we will adopt the same approach 
here in estimating the contribution of the four fermion sources to the neutron EDM. 

To proceed, it is convenient to first introduce a redefined set of Wilson coefficients, so 
that 

6C = Cbd{'m'b)'m'bbbdi'yrjd + Cdh{'rni))mi,ddhi'~^^h, (A. 65) 

at the 6-quark mass scale, where C = C /mi,. Below this scale, the 6-quarks can be integrated 
out leading to dimension-7 quark-gluon operators. Making use of the standard low energy 
theorems we obtain 

^ _?^A!^m!^di^^dGG - ^^f^p!^ddGG. (A.66) 

The evolution of the coefficients to the hadronic scale ~ 1 GeV requires knowledge of 
the anomalous dimensions of these dimension-7 operators. Since (with hindsight), the 
overall contribution of these operators is small, we will approximate this by the anomalous 
dimension of qq, recalling that agGG is RG-invariant. Thus, we run and Cdb to // ~ 
1 GeV, using -fasqqGG ~ 4/9. 

To gain some intuition about the dependence of d„ on C(//), we first utilize "naive 
dimensional analysis" [54, 55]. The quark-gluon operators are of dimension-7, so the reduced 
coupling is M^/(47r)^ with M ~ 47r/^. The corresponding estimate for dn is given by e/M 
multiplied by the reduced coupling of the source, and so we find 

]d„(NDA)| ~ e 2.1 X 10"=^ GeV' (^^(/x) + 0.75 CM) , (A.67) 

on using the usual matching condition ag = 47r, i.e. at a very low (apparent) scale. As it 
turns out, our final estimate is quite close to this, assuming we choose the matching scale 
not too far from 1 GeV. 

We now turn to a sum-rule analysis. As alluded to above, a direct calculation appears 
intractable due to the presence of unknown condensates. Therefore, we will pursue instead 
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Figure 7: The leading contribution to the 75 structure induced by the CP-odd source. 

an indirect method first considered by Bigi and Uraltsev [76], and used recently within 
QCD sum-rules to estimate dn induced by the Weinberg operator GGG. One considers the 
75-rotation of the nucleon wavefunction induced by SC, leading to an estimate for dn in 
terms of the corresponding rotation of the anomalous magnetic moment 

{N\6C\N) 

In the present context, we require the ratio of the leading contributions to the structures 
1 and in the mass sum-rule correlator. We have 

I d'xe'P-^{rin(0)fin(x)) = Y^(^g) [p'ln( V)l + fsip')^^] , (A.69) 

where rjn is the neutron current, which for the mass sum-rule we take at the "lofi^e point" 
[77] , fs is the contribution induced by the CP-odd souce 6jC that is to be determined, and 
we have retained only the leading order contribution to the structure 1. 

Fortunately, in this case the leading contribution enters in a tractable form without the 
need for additional perturbative vertices - i.e. the gluon structure of the source SC can 
be extracted at the nonperturbative level. The relevant class of diagrams is exhibited in 
Fig. 7. Using standard methods, we calculate the quark propagators in the background of 
the two sources in 5C making use of the known guark-gluon condensate^: 

{qgs{Ga)q) = ml\{qq)\ ~ 0.8GeY'{qq). (A.70) 

One then obtains, 

Mp') = (^"'^ + ^-^^Cdb) , (A.71) 

and we proceed by performing a Borel transform and comparing the two terms in (A.69) 
at the neutron mass scale in terms of the Borel parameter M. This fixes the ratio of the 

®We have introduced a tilde atop mo in order to distinguish this QCD condensate from the mass of the 
SUSY scalars. 
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two structures in (A. 69) which is the quantity that enters the estimate (A. 68), leading to 
the result, 

\dn\ ~ \l^n\-^,{qq) (Cm + 0.75^^^) - e 2.6 x IQ-^GeV^ [Cm + 0.75Cab) , (A.72) 

which necessarily comes with a precision estimate not better than C(100%). On re- 
expresing this result in terms of Cm and Cdb, we obtain the result quoted in (4.57). 

Finally, for completeness, we note that the operators Cdd, Cds, Cgd and Css can be 
generated via H — A mixing. Their contribution to the neutron EDM can be estimated as 
[9, 35]: 

^ dJ^ ^ 2 X 10-^GeV, ^^^^ 2 x IQ-^GeV. (A.73) 

Cdd Cgd Cds Css 



B Mercury EDM induced by four-fermion operators 

The main source for the Schiff moment of the mercury nucleus arises through a pion- 
exchange contribution with CP violation in the pion-nucleon vertex, 

C.NN = tLNr^^Nn- + £LnNt:'' + ^S^(7VT«A^7r" - ZNt^'Nt:''), (B.74) 

of which the contribution of ^^jvAf the most important [67, 66]. We recall here that 
the main technical difficulty in the derivation of this couphng, as induced by color EDM 
operators, is the existence of two compensating diagrams. The direct chiral commutator 
[Jo5, g(Gcr)75g] tends to cancel against a pion re-scattering diagram (see Ref. [67] for de- 
tails). The resulting expression for (?^^jvAf given by the following linear combination of 
nucleon matrix elements for qgs{Ga)q and qq: 

du dd/j^ji- \ 7 \, ~ 2/- 
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{N\ugs{Ga)u + dgs{Ga)d - mi{uu + dd)\N). (B.75) 



This vanishes if {N\qgs{Ga)q\N) — m^{N\qq\N) , a relation that holds in the vacuum 
factorization approximation. Thus, in order to obtain a non-zero result, one has to go 
beyond the factorization approximation. A full QCD sum-rule analysis [67] produces a 
preferred range for this coupling with the "best" value given by 

-g^^l^ = 2 X 10^^ (J„ - dd) cm-\ (B.76) 

Combined with the relevant nuclear, S{g^%j^), and atomic, diig{S), pieces of the calculation, 
this expression enters the first line of Eq. (4.61). 

After this brief review of the color EDM contributions to g^pfj^f, for large tan/3 we also 
need to generalize this derivation to include gi%N{Cqiq2)- Both diagrams, the direct chiral 
commutator and the pion re-scattering diagram, are present in this case as well. Luckily, 
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in contrast with the case of CEDM sources, vacuum factorization does now produce a non- 
zero result. Combining these two diagrams, we find a contribution that corresponds to the 
direct factorization of the pseudoscalar bihnear qi^^q via the chiral commutator with the 
pion. The final expression for g^)fpf takes the form, 

9^L = ^{N\CMdd + C,ass + C,M\N) 

= -8xl0-2GeV2f^^:^ + 3.3/^^ + ^(l-0.25«)V (B.77) 

V md rus rrib J 

Note that within the factorization approximation Cat, Cds, Cgs and Cbi, do not contribute 
to ^^TViv- The final expression (B.77) depends on the same matrix element of ss over the 
nucleon as previously encountered in the calculation of C5, and parametrized by k. 
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